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We investigate the magnetic properties of three Mng single molecule magnets by means of inelastic 
neutron scattering and frequency domain magnetic resonance spectroscopy. The experimental data 
reveal that small structural distortions of the molecular geometry produce a significant effect on the 
energy level diagram and therefore on the magnetic properties of the molecule. We show that the 
giant spin model completely fails to describe the spin level structure of the ground spin multiplets. 
We analyze theoretically the spin Hamiltonian for the low spin Mn6 molecule (S — 4) and we 
show that the excited S multiplets play a key role in determining the effective energy barrier for the 
magnetization reversal, in analogy to what was previously found for the two high spin Mn6 (S = 12) 
molecules [S. Carretta et ai, Phys. Rev. Lett. 100, 157203 (2008)]. 



PACS numbers: 75.50.Xx, 78.70.Nx, 33.35.+r, 75.60.Jk 
I. INTRODUCTION 



Single Molecule Magnets (SMMs) have been the 
subject of intense research activity since the first and 
mostly studied one, Mni 2 -ac, was reported 1 -. These 
metal-organic clusters are usually characterized by a 
large spin ground state S and an easy-axis anisotropy 
which determines the Zero-Field Splitting (ZFS) of the S 
state sublevels. The resulting magnetic bistability makes 
them interesting for magnetic storage applications due 
to their potential to shrink the magnetic bit down to the 
size of one single molecule. Until recently and despite 
the common efforts of chemists and physicists to find 
suitable systems that could retain the magnetization for 
a long time at non cryogenic temperatures, Mni2-ac was 
the system showing the 'highest' blocking temperature 
(3.5 K) and anisotropy barrier (74.4 K)2. The relaxation 
time in the classical regime follows the Arrhenius law: 
t = T exp(U/kBT) (Ref. H). According to this, there 
are two key points that have to be considered for the 
realization of an ideal SMM. First of all, the anisotropy 
barrier, given to a first approximation by U ~ \D\S (D 
is the axial anisotropy parameter), has to be sufficiently 
high. This is to prevent the reversal of the magnetization 
via a classical thermally activated multistep Orbach 
process mediated by spin-phonon interactions. This 
can be achieved by the simultaneous increase of D and 
S, two variables that are intrinsically linked together 4 . 
Secondly, the pre-exponential factor tq in the Arrhenius 



law has to be large. This factor is dominated by the 
time necessary to climb the upper states in the energy 
level diagram, and is proportional to D~ 3 (Ref. y). 
In addition to the classical relaxation mechanism, the 
quantum tunneling of the magnetization (QTM) that 
characterizes the spin dynamics of SMMs, has to be 
taken into consideration and minimized for magnetic 
data storage application, since it provides a shortcut for 
the relaxation of the magnetization. 
Therefore, to engineer SMMs able to retain the mag- 
netization for long time it is crucial to control all the 
different mechanisms that provide a relaxation path for 
the system. Recently we succeeded in the synthesis of 
a new class of Mn 3+ -based clusters that contributed 
in raising the anisotropy barrier and has served as a 
good model system to study the factors involved in the 
relaxation mechanism 5 ^. 

This class consists of hexanuclear Mn 3+ clusters (from 
now on Mn@) which, despite the generally similar 
nuclear structure, display a rich variety of spin ground 
states and anisotropy energy barriers^—. The six 
Mn 3+ ions are arranged in two triangles, with dominant 
ferromagnetic (FM) exchange interaction between the 
two triangles and FM or antiferromagnetic (AFM) 
interactions within the two triangles. It has been found 
that the nature of the intra-triangle exchange interaction 
can be switched from AFM to FM by substituting the 
organic ligands bridging the Mn 3+ ions, leading to a 
change of the ground state from a low spin (S = 4) 
to a high spin (S — 12) 7 . Furthermore, deliberately 
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targeted structural distortions have been successfully 
used to tune the values of the exchange interactions^. 
The isotropic exchange interactions, and consequently 
the overall anisotropy barrier—, is thus found to be very 
sensitive to the structural details. This has been also 
demonstrated using an alternative method for distorting 
the molecule, that is by applying external hydrostatic 
pressure and correlating the structural changes with the 
magnetic behavior—. It is therefore quite important 
to determine the exchange interactions for different 
structures to deduce magneto-structural correlations. 
This information can be then used to engineer new 
clusters with selectively modified molecular structures 
that match the optimized conditions for the desired 
magnetic properties. 

We have investigated three members of the 
family of Mng clusters, with chemical formu- 
las [Mn 6 2 (sao) 6 (0 2 CMe)2(EtOH)4]-4EtOH (1), 
[Mn 6 2 (Et-sao) 6 (0 2 CPh) 2 (EtOH) 4 (H 2 0) 2 ]-2EtOH (2) 
and [Mn 6 2 (Et-sao) 6 (0 2 CPh(Me) 2 ) 2 (EtOH) 6 ] (3)£&M. 
All molecules display very similar structures consisting 
of six Mn 3+ ions (s = 2) arranged in two staggered 
triangular units (see Fig. [1} related by an inversion 
centre. 




FIG. 1: (color online) Core of molecules (1) (left) and (3) 
(right) showing at the bottom the difference in torsion angles 
(ai, Qf2 and 03). Color scheme: Mn, orange, O, red, N, blue. 
H and C ions are omitted for clarity. 

The only major structural difference between the three 
clusters resides is the steric effect of the organic ligands 
used in proximity to the transition metal ions. How- 
ever, despite having very similar structures, the three 
molecules have very different magnetic properties. The 
coupling between the magnetic ions occurs via superex- 
change pathways involving oxygen and nitrogen ions and 
is found to be extremely sensitive to intramolecular bond 
angles and distances. The particular arrangement of the 
magnetic ions provides exchange couplings lying in the 
cross-over region between AFM and FM. For this reason, 
even small structural distortions have tremendous impact 
on the magnetic properties of the system. For example, 
while the coupling between the two triangles is ferro- 
magnetic for all molecules, the intra-triangular coupling 
changes from antiferromagnetic in (1) to ferromagnetic 



in (2) and (3) due to a 'twisting' of the oximate linkage. 
This results in a 'switching' of the total spin ground state 
from S' = 4to5 = 12. Systematic synthesis and stud- 
ies of various members of the Mng family have revealed 
that the nature of the coupling is extremely sensitive to 
the intra-triangular Mn-O-N-Mn torsion angle s 10 ' 11 (see 
Fig. [TJ). There is a critical value for the torsion angle of 
30.85 ± 0.45°, above which the pairwise exchange inter- 
action switches from antiferromagnetic to ferromagnetic, 
while a further enhancement of the angle increases the 
strength of the FM interaction. This effect has been in- 
terpreted in terms of the particular arrangement of the 
manganese d 2 2 orbitals with respect to the p-orbitals of 
the nitrogen and oxygen ions. A large (small) Mn-O- 
N-Mn torsion angle results in a small (large) overlap 
between the magnetic orbitals giving rise to ferromag- 
netic (antiferromagnetic or weak ferromagnetic) superex- 
change interactions^. 

Molecules (2) and (3) have the same spin ground state 
S = 12 but very different effective energy barriers (U e g « 
53 K for (2) and U cS « 86.4 K for (3)). This differ- 
ence was found to be closely related to the exchange 
interactions^. 

In order to understand this rich variety of behaviors, 
we performed a detailed spectroscopic characterization 
of the three molecules using inelastic neutron scatter- 
ing (INS) and frequency domain magnetic resonance 
(FDMR). FDMR is only sensitive to transitions with 
a predominantly intra-multiplet character, according to 
the selection rules AS — 0,AMg = ±1. In contrast, 
in INS both inter- and intra-multiplet transitions can be 
observed (AS = 0,±1,AM S = 0,±1). Thus, the com- 
bination of the two techniques allows assignment of all 
observed excitation s 16 ' 17 . 

The determination of the model spin Hamiltonian pa- 
rameters enabled us to estimate theoretically the effec- 
tive energy barrier for the low spin molecule (1). Sim- 
ilarly to what we previously reported for the two high 
spin molecules (2) and (3), the results on (1) show how 
the presence of low-lying excited spin multiplets plays a 
crucial role in determining the relaxation of the magne- 
tization. 

In conventional systems, the effects of S-mixing can be 
effectively modeled by the inclusion of fourth order zero- 
field splitting parameters in the giant spin Hamiltonian^. 
Here we will show that this Hamiltonian is completely 
inadequate for the description of the spin state energy 
level structure. 



II. EXPERIMENTAL METHODS 

Non-deuterated polycrystalline samples were synthe- 
sized according to published methods^. 

FDMR spectra were recorded on a previously described 
quasi-optical spectrometer^ which employs backward 
wave oscillators as monochromatic coherent radiation 
sources and a Golay cell as detector. Sample (1) proved 
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to deteriorate rapidly upon pressing and over time. 
Therefore, the FDMR measurements on (1) were per- 
formed on loose microcrystalline material (348 mg) held 
between two quartz plates. In this unconventional mea- 
surement, the detector signal was recorded as function of 
frequency at different temperatures. Extreme care had 
to be taken to prevent the slightest positional changes 
of sample and equipment, which changes the standing 
wave pattern in the beam, precluding normalization. The 
normalized transmission was calculated by dividing the 
signal intensity at a given temperature by that at the 
highest temperature (70K). Sample (2) and (3) deterio- 
rate to a lesser extent and FDMR spectra were recorded 
on pressed powder pellets made by pressing ca. 250 mg 
of the unground sample, with ca. 50 mg n-eicosane (to 
improve pellet quality) into a pellet. All spectra were 
simulated using previously described software.— 

INS experiments were performed using the multi disc- 
chopper time-of- fight spectrometers V3/NEAT at the 
Helmholtz-Zentrum Berlin fur Materialien und Energie 
(HZB, Berlin, Germany) and IN5 and IN6 at the Institute 
Laue Langevin (Grenoble, France). The samples were in- 
serted into hollow cylindric shaped Aluminum containers 
and mounted inside a standard orange cryostat to achieve 
a base temperature of 2 K. A vanadium standard was 
used for the detector normalization and empty can mea- 
surements were used for the background subtraction. 

III. THEORETICAL MODELING AND 
EXPERIMENTAL RESULTS 

The experimental data have been modeled using both 
the giant spin Hamiltonian (GSH), which considers the 
ZFS of the ground state multiplet only, and the micro- 
scopic spin Hamiltonian, which treats isotropic exchange 
and single-ion ZFS at the same level. Including only 
ZFS terms, the giant spin Hamiltonian for a spin state S 
reads: 

H s = D s S 2 z +E s (S 2 x -S 2 v ) + B° 4 dl (1) 

where D$ and E$ are second order axial and transverse 
anisotropy, respectively, and B° is the fourth order axial 
anisotropy, with O4 the corresponding Stevens operator. 
The microscopic spin Hamiltonian includes an isotropic 
exchange term for each pairwise interaction and single 
ion ZFS terms for each ion: 

+Ci(25-30a(s + l))s2(*)). (2) 

where s(i) are spin operators of the i th Mn ion. The first 
term is the isotropic exchange interaction, while the sec- 
ond and third terms are the second and fourth order axial 
single-ion zero-field splitting, respectively (the z axis is 



assumed perpendicular to the plane of the triangle). 
The spin Hamiltonians have been numerically diagonal- 
ized by exploiting the conservation of the z-component of 
the molecular total spin and the exchange and anisotropy 
parameters have been varied to obtain a best fit of the 
experimental data. 

A. Mn 6 (1) (S = 4) U ea w 28 K 

Sample (1) was the first reported member of the Mng 
family 14 . The building block of the molecule is the [Mng + 
O] triangular unit where Mn2 pairs, bridged by the NO 
oximc, form a -Mn-O-N-Mn- moiety (Fig. [2]). 




FIG. 2: (color online) Structure of the Mns (1) molecular 
core. The Mn 3+ ions are located at the vertices of two oxo- 
centered triangles. Ions Mnl, Mn2, Mnl' and Mn2' are in 
octahedral geometry and ions Mn3 and Mn3' in square pyra- 
midal geometry, as highlighted in filled and striped orange 
(left figure). Color scheme: Mn, orange, O, red, N, blue. H 
and C ions are omitted for clarity. 

The Mn-O-N-Mn torsion angles within each triangle 
are 10.7°, 16.48° and 22.8°, giving rise to a dominant 
antifcrromagnetic exchange coupling 11 -. The two trian- 
gular units are coupled ferromagnetically, resulting in a 
total spin ground state of S = 4. Four out of the six 
metal ions (Mnl, Mn2, Mnl', Mn2') are six-coordinate 
and in distorted octahedral geometry (MnOsN), with the 
Jahn- Teller axis almost perpendicular to the plane of the 
triangle, while the two remaining ions (Mn3, Mn3') are 
five-coordinate and in square pyramidal geometry (see 
Fig. [2]). The effective energy barrier was determined 
from AC susceptibility measurements to be U c s = 28 K, 
with To = 3.6 x 10~ 8 s (Ref. 14). From the effective en- 
ergy barrier an estimate of D w —0.15 meV was derived. 

We performed INS and FDMR measurements to char- 
acterize the ground multiplet and to identify the posi- 
tion of the lowest-lying excited states from which we de- 
termine the effective exchange interaction and the zero- 
field splitting parameters. Figure [3] shows the FDMR 
spectra recorded on 350 mg unpressed powder of (1). 
The most pronounced feature is the resonance line at 
1.803(7) meV, while much weaker features can be ob- 
served at 1.328(1) meV and 1.07(1) meV. The intensity of 
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FIG. 3: (Color online) (a) FDMR spectra of unpressed poly- 
crystalline powder of (1) recorded at various temperatures. 
The intensity of the higher-frequency resonance line decreases 
with temperature, while that of the lower-frequency lines in- 
creases up to 30 K, beyond which it decreases again. Dotted 
lines indicate resonance lines due to impurities, (b) Expanded 
view of the low frequency part of the 30 K spectrum, (c) Ex- 
perimental and fitted spectrum at T — 30 K using the GSH 
approximation. Note the logarithmic scale in (a) and (c). 



the higher-frequency line is strongest at lowest tempera- 
ture, proving that the corresponding transition originates 
from the ground spin multiplet. The lower-frequency 
lines have maximum intensity at around 30 K. No fur- 
ther features were observed between 0.5 and 3 meV. The 
intense resonance line shows two shoulders to lower ener- 
gies, which are much stronger in pressed powder samples 
and also increase with the age of the sample. This be- 
havior is mirrored by the development of a pronounced 
asymmetric lineshape in INS studies on older samples. 
We attribute these shoulders to microcrystalline parti- 
cles that have suffered loss of lattice solvent, which leads 
to small conformational changes and this alters the ZFS 
and exchange parameters. We discard the possibility of 
isomers with different orientations of the Jahn- Teller dis- 
tortion axes, as observed for Mnisri, because we see no 
signature of different isomers in the AC susceptibility. 
We also discount the possibility of closely spaced tran- 
sitions, as observed in the Fei3 cluster 22 , because the 



intratriangle exchange interactions are not equal. 

The higher frequency resonance line is attributed 
to the transition from the \S = 4, Ms = ±4) to 
\S = 4, Ms ± 3) states. INS measurements have found 
to be necessary to unambiguously identify the origin of 
the lower frequency transitions (see below). Assuming 
that these transitions are transitions within the ground 
multiplet, a fit of the giant spin Hamiltonian ZFS pa- 
rameters (Eq. [T]) to the observed resonance line energies 
yields D s=4 = -2.12 ± 0.03 cm" 1 (-0.263 ± 0.004 meV) 
and Bl = +(1.5 ± 0.5) x 10~ 4 cm" 1 (1.24 ± 0.06 x 10~ 5 
meV). This ground state D$- value is much larger than re- 
ported spectroscopically determined Ds-paramcters for 
other manganese SMMs, e.g. Ds=io = —0.457 cm -1 
for Mni 2 Ac 23 , D s=1 t_ = -0.247 cm" 1 for Mn^, or 

Ds=e = —1.16 cm -1 for Mn^Zna 24 . The main rea- 
son for this large D- value is the fact that the projec- 
tion coefficients for the single ion ZFS onto the clus- 
ter ZFS are larger for spin states with lower S (Ref. 
|25| ). The determined Ds=a = —2.12 cm -1 value for 
(1) is in excellent agreement with that found from DFT 
calculations (D = —2.15 cm -1 ) 26 . The expected en- 
ergy barrier toward relaxation of the magnetization cal- 
culated from the found spin Hamiltonian parameters is 
t^thcor = \D\S 2 = 48.8 K, which is much larger than 
the experimentally found U e s ~ 28 K, indicating that 
more complex relaxation dynamics characterize this sys- 
tem, in analogy to what has been found for the Mn 6 
S = 12 compounds^. The linewidth of the 1.33 meV line 
is slightly larger than that of the 1.80 meV line (48 /ieV 
versus 41 /ieV), which can indicate the presence of more 
than one excitation. The simulated spectrum agrees very 
well for the higher-frequency resonance line (note that 
the intensity is not rescaled), while the lower- frequency 
line is much weaker in the experiment than from the fit. 
This can be tentatively attributed to the presence of low- 
lying excited states as observed previously for Mn^Ac 27 . 
To determine the energy of excited spin states and iden- 
tify the origin of the low frequency resonances we resorted 
to INS, the technique of choice to directly access inter- 
multiplet excitations. 

The INS experiments were performed on w 4 g of non- 
deuterated polycrystalline powder of (1), which was syn- 
thesized as described in Refll4 For our measurements 
we used incident neutron wavelengths ranging from 3.0 
A to 5.92 A with energy resolution between 50 /xeV and 
360 fxeV. 

Figure @] (a) shows the INS spectra for an incident wave- 
length of 4.6 A collected on NEAT (210 fi eV full width at 
half maximum (FWHM) resolution at the elastic peak) . 
At T = 2 K, only the ground state is populated and 
therefore all excitations arise from the ground state dou- 
blet \S = 4, Ms = ±4). We observed a strong transition 
at 1.77(2) meV, which we assign to the intra-multiplet 
transition to the \S — AMs — ±3) level, in agreement 
with FDMR results (see above). One further excitation 
was observed at higher energy at 2.53(1) meV. 

At T — 20 K, we detected additional excitations 
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at 1.05(1) meV and 1.31(1) meV, which must be due 
to transitions from excited states. All peaks in the 
INS spectra show a very unusual asymmetric line-shape, 
which we assign to lattice solvent loss (see above). 
From the comparison of INS data with the FDMR re- 
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FIG. 4: (Color online) (a) INS spectra of (1) with an incident 
wavelength of A = 4.6 A (NEAT) for T = 2 K (blue circles) 
and T = 20 K (red squares). The continuous lines represent 
the spectra calculated assuming a dimer model for the spin 
Hamiltonian (Eq. 3). (b) Q-dependence of first intra- (green 
circles) and inter-multiplet (black squares) transitions mea- 
sured on IN6 for A = 4.1 A and T=2 K. Continuous lines 
represent the calculated Q-dependence using the dimer spin 
Hamiltonian Eq. @ (assuming a dimer distance of R = 5.17 
A, which corresponds to the distance between the centre of 
the two triangles). 

suits, we can deduce that the excitation at 2.53 meV 
has a pure inter-multiplet origin, being absent in the 
FDMR spectra (see Fig. [3]). This is also confirmed 
by the Q-dependence of the scattering intensity of the 
observed excitations. Figure 0] (b) shows this depen- 
dence for the \S = 4, M s = ±4} \S = 4 M s = ±3) and 
\S = 4, M s = ±4) -4 \S = 3, M s = ±3) transitions. A 
characteristic oscillatory behavior has been observed for 
the Q dependence of the inter-multiplet INS transition 
(black squares), which presents a maximum of intensity 
at a finite Q value (that is related to the geometry of 
the molecule) , and decreasing intensity as Q goes toward 
zero. This Q dependence is typical for magnetic clusters 



and reflects the multi-spin nature of the spin state o 28 ' 29 . 
By contrast, the intra-multiplet excitation (green circles) 
has maximum intensity at Q = 0, as expected for a tran- 
sition with AS* = 0, and the intensity decreases with 
increasing Q, following the magnetic form factor. 

The INS data directly reveal the presence of low ly- 
ing excited multiplets. Indeed, the difference in energy 
between the lowest and the highest energy levels of the 
anisotropy split 5 = 4 ground state is given, as a first 
approximation, by \D\S 2 =4.2 meV. The presence of an 
inter-multiplet excitation at only 2.53 meV energy trans- 
fer, therefore below 4.2 meV, indicates that the first ex- 
cited 5 multiplet lies within the energy interval of the 
anisotropy split 5 = 4 state. This suggests that the ob- 
served low energy excitations are possibly not pure intra- 
multiplet transitions, but are expected to originate from 
the 5 = 4 ground state and from the first excited 5 multi- 
plet. Therefore the exact assignment of those excitations 
requires a more accurate analysis beyond the GSH ap- 
proximation. Indeed, one fundamental requirement for 
the validity of the GSH approximation, i.e. an isolated 
ground state well separated from the the excited states, 
is not fulfilled and 5 is not a good quantum number to 
describe the ground state of the molecule. To model the 
data it is thus necessary to use the full microscopic spin 
Hamiltonian of Eq. [2j 

Given the low symmetry of the triangular units in (1), 
the number of free parameters in Eq. [2] would be too 
large to obtain unambiguous results, considering the low 
number of experimentally observed excitations. Hence, 
we have chosen to describe the low-energy physics of (1) 
by a simplified dimer model, an approximation which has 
already previously been adopted for (2) and (3) (see Ref. 
l3fjl ). More specifically, the two triangular units are de- 
scribed as two ferromagnetically-coupled 5 = 2 spins, 
which also experience an effective uniaxial crystal-field 
potential: 

H dimer = J(S A ■ 8 B ) + d(S% z + S 2 B z ). (3) 

The spin Hamiltonian has been diagonalized numerically 
and the J and d parameters have been varied to obtain 
a best fit of the experimental data. The position of the 
peak at 1.77 meV does not depend on the exchange in- 
teraction, therefore its position sets the value of the axial 
anisotropy d parameter. Given the d parameter, a fit of 
the position of the peak at 2.53 meV sets the isotropic 
exchange parameter J. 

The best fit of the experimental data is obtained with 
J = —0.19 meV and d = —0.59 meV. The calculated en- 
ergy level scheme is reported in Fig. [5] (left), where the 
comparison with the energy level diagram in the GSH 
approximation is also reported (right). The value of S e ff 
(where (5 2 ) := 5 e g(5 e ff+l)) is labeled in color and shows 
that the first 5 = 3 excited state is completely nested 
within the 5 = 4 ground state. From Fig. [5] it is also 
clear that the GSH model does not account for a number 
of spin states different from the ground state 5 = 4 mul- 
tiplet at low energy. Furthermore, the assignment of the 
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groups on the carboxylate ligand has increased the non- 
planarity of the Mn— N— O— Mn moieties further, giving 
torsion angles of 39.08°, 43.07° and 34.88° 8 . The re- 
sult is that the weakest ferromagnetic coupling is signif- 
icantly stronger for (3) compared to (2). Using a single 
J model (e.g. assuming that the intra- and inter-triangle 
exchange couplings are equal), Milios et al. fitted the 
DC susceptibility data for molecules (2) and (3) and ob- 
tained: J(2)= -0.230 meV and J(3)= -0.404 meV, 
respectivel y 10 ' 12 (in our notation for the spin Hamilto- 
nian). 

In spite of the fact that both (2) and (3) have S = 12 



FIG. 5: (Color online) Calculated energy level diagram for 
molecule (1). The level scheme on the left side is calculated 
using the dimer model spin Hamiltonian in Eq. [3] The color 
maps S e e , where (S 2 ) := S e e(S e s+ 1). The black dashed line 
corresponds to the observed value of U e g = 28 K. The black 
arrows indicate transitions which contribute to the observed 
peak in the INS and FDMR spectra at E ~ 1.33 meV (see 
text for details). The level diagram on the right has been 
calculated using the GSH approximation (Eq. [TJ. 



observed excitations can be misleading if considering the 
GSH approximation only. For example, using the GSH 
model, the observed peak at 1.33 meV can only be at- 
tributed to a pure intra- multiplet excitation from 4, ±3) 
to 4, ±2), whilst using Eq. [3J it is found to be a su- 
perposition of several inter-multiplet and intra-multiplet 
transitions (indicated by arrows in Fig. [SJ. The GSH 
approximation fails to describe the low energy level dia- 
gram of the molecule and consequently fails to describe 
the relaxation of the magnetization. Indeed, the presence 
of excited states nested within the ground state multiplet 
has a significant effect on the relaxation dynamics, as dis- 
cussed in section HVl 



B. Mn 6 (2) U eft » 53 K vs. Mn 6 (3) U eS ~ 86.4 K 

Introducing sterically more demanding oximate ligands 
results in a twisting of the Mn-N-O-Mn torsion angle 7 , 
which causes switching of the intra-triangle exchange in- 
teractions from antiferromagnetic to ferromagnetic, re- 
sulting in a large increase of the spin of the ground state 
from S = 4 to S = 12. Here, we study two ((2) and 
(3), respectively) of the many published derivatives of 
these S — 12 Mn6 clusters. Compound (2) has un- 
dergone two structural changes compared to (1). First 
of all, the distance between the phenolato oxygen and 
the two square pyramidal Mn 3+ ions has decreased from 
rj 3.5 A to « 2.5 A, thus all Mn 3+ ions arc now m six- 
coordinated distorted octahedral geometry (see Fig. [5]). 
Secondly, the torsion angles of the Mn— N— O— Mn moi- 
eties has increased strongly with respect to those in (1), 
being 38.20°, 39.9° and 31.26°, compared to 10.7°, 16.48° 
and 22.8° for (1). In (3), the introduction of two methyl 
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FIG. 6: (Color online) Structure of the Mn6 (2) molecu- 
lar core. The Mn 3+ ions are located at the vertices of two 
oxo-centered triangles. All Mn ions are in octahedral geome- 
try and the octahedra are highlighted in orange (left figure). 
Color scheme: Mn, orange, O, red, N, light blue. H and C ions 
are omitted for clarity. On the right, a schematic represen- 
tation is given, together with the exchange coupling scheme 
adopted for the spin Hamiltonian calculations. 

ground states and similar geometrical structures, radi- 
cally different effective energy barriers towards the relax- 
ation of the magnetization were observed, being U c g ~ 53 
K for (2) and U eS m 86.4 K for (3). Here, we aim to 
understand this difference by an in-depth study of the 
energy level structure by means of FDMR and INS. 
Figure [7] shows FDMR spectra recorded on a pressed 
powder pellet of (2) at different temperatures. The 
baseline shows a pronounced oscillation, which is due 
to Fabry-Perot-like interference within the plane-parallel 
pellet^. The oscillation period and downward slope to 
higher frequencies are determined by the thickness of the 
pellet and the complex dielectric permittivity, which were 
determined to be e' = 3.01 and e" = 0.049, values typical 
for molecular magnet samples. In addition, five reso- 
nance lines are observed which we attribute to resonance 
transitions within the S — 12 multiplet. Thus, the high- 
est frequency line is assigned to the 1 12, ±12) — >• 1 12, ±11) 
transition, and so on. The lines are much narrower 
(11 /ueV FWHM) than those observed for other SMMs, 
e.g. 23 /ieV FWHM for Mn^Ac. The fit procedure 
showed that the lines are inhomogeneously broadened 
and best described by Gaussian lineshapes. The small 
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Energy (meV) 
0.75 1.00 



TABLE I: INS and FDMR peak positions of the observed 
excitations for (2) and (3) (in meV). 
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FIG. 7: (Color online) (a) FDMR spectra recorded on a 
pressed powder pellet of 2 at different temperatures. At the 
lowest temperature, the highest-frequency line has highest in- 
tensity. The other transitions are indicated by vertical lines, 
(b) 10K FDMR spectrum (symbols) and best fit using the 
GSH (Eq. [[]), with D = -0.368 cm" 1 , and Bt = -4.0 x 10" 6 



linewidth indicates that distributions in ZFS parameters 
(D-strain) are small in these samples. A fit of the GSH 
parameters (Eq. [1]) to the observed resonance frequen- 
cies, yields D s=12 = -0.368 cm" 1 (0.0456 meV) and 
Bl = -4.0 x 10" 6 cm" 1 (4.96xl0" 7 meV) best parame- 
ter values. The theoretical energy barrier calculated from 
these ZFS parameters is t/thcor = 76 K, which is much 
larger than the experimentally found U c s ~ 53 K, indi- 
cating that the molecule can shortcut the barrier in some 
way. The ZFS values are in themselves not remarkable, 
and close to those reported for other manganese clusters 
with similar ground state spins, e.g. Ds=w = —0.457 
cm" 1 for Mn 12 A(M, D s= it_ = -0.247 cm" 1 for Mnc^l 
Interestingly, the fourth order axial ZFS is an order of 
magnitude smaller than for Mn^Ac. This type of ZFS 
is currently accepted to parametrize effects of mixing be- 
tween spin multiplets (S'-mixing)^, which would mean 
that S'-mixing is only limited, contrary to expectation. 
However, the fit does not simulate the resonance line po- 
sitions and intensities satisfactorily, which is in contrast 
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to the situation for other molecular nanomagnets that 
feature strong S'-mixing, e.g. N1 4 16 ' 32 . Therefore, the in- 
vestigated Mn 6 SMM represents an example where the 
giant spin model cannot satisfactorily describe FDMR 
spectra, and it will be shown below that this is due to a 
complete breakdown of the giant spin model. It will also 
be shown that the resonance line at 0.80 meV is due to a 
transition within the 5=11 excited multiplet. However, 
removal of this resonance line does not result in a better 
fit. The calculated line intensities are much larger than 
those experimentally found, especially for the highest- 
frequency lines. This we attribute to a combination of 
parasitic radiation in the cryostat, and the presence of 
many more states than taken into account by the giant 
spin model, which decreases the relative population for 
any given state. 

Similar FDMR results were obtained for (3) (Fig. [5]), 
and six sharp resonance lines were observed. A fit of 
the GSH parameters to the observed resonance line po- 
sitions yields the following values: D = —0.362 ± 0.001 
cm" 1 (-0.0449 meV), and Bl = -6.0 ± 0.4 x 10" 6 cm" 1 
(-7.4 xlO -7 meV). The simulated spectrum matches the 
experiment much more closely for (3), especially for the 
high-frequency lines. Interestingly, the theoretical energy 
barrier (L/th oor = 75 K) is virtually the same as for (2), 
but smaller than the experimentally found energy barrier 
(U e s = 86 K). This unprecedented finding means that the 
magnetization relaxation must involve states that do not 
belong to the ground spin multiple^. Again, we turn to 
INS to determine the positions of the excited spin multi- 
plets, which will allow full characterization of the system. 
Figures (|TJ£i) and (|TJp) show the high resolution INS ex- 
perimental data for compounds (2) and (3), respectively, 
collected on IN5 with an incident wavelength of 6.7 A (53 
fieV FWHM resolution at the elastic peak) . At the lowest 
temperature T — 2 K only the ground state is populated 
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FIG. 8: (Color online) (a) FDMR spectra of unpressed poly- 
crystalline powder of (3) recorded at various temperatures, 
(b) Experimental and fitted spectrum at T = 10 K. 



and, due to the INS selection rules, only transitions with 
AS = 0, ±1 and AM = 0, ±1 can be detected. The low- 
est energy excitation can be thus easily attributed to the 
intra-multiplet transition from the \S = 12, Ms — ±12} 
ground state to the \S = 12, Ms = ±11) first excited 
level. The position of this intra-multiplet excitation is 
found to be at about the same energy in both compounds, 
i.e. ~ 1.1 meV, indicating only small differences in the 
anisotropy of the two systems. In contrast the first inter- 
multiplet S — 12 — > S = 11 excitation at about 1.41 
meV in compound (2) is not visible in the spectra at 6.7 
A of compound (3). This can be understood looking at 
the data at higher energy transfer, collected with an in- 
cident wavelength of 3.4 A (see Figs. (fTUk ) and (fTOb)). 
Indeed the first inter-multiplet excitation is considerably 
raised in energy in compound (3) with respect to com- 
pound (2), from 1.41 meV to 1.87 meV. This gives a 
direct evidence of an increase of the isotropic exchange 
parameters, while the anisotropic parameters are approx- 
imately the same for both molecules. The INS spectra 
collected at a base temperature of 2 K, enabled us to 
directly access the whole set of intra-multiplet and inter- 
multiplet transitions allowed by the INS selection rules 
in both compounds. By raising the temperature to 16 




0.5 1 
Energy transfer (meV) 



FIG. 9: (Color online) INS spectra collected on IN5 with 
incident wavelength of 6.7 A at T = 2 K (blue circles) and 
T = 12 K (red squares), (a) for sample (2) and (b) for sample 
(3). The spectra calculated with the parameters listed in 
Table ITT1 are shown as continuous lines. 



K the intensity of the magnetic peaks decreases, thus 
confirming their magnetic origin. A total of five inter- 
multiplet excitations for compound (2) toward different 
5 = 11 excited states can be detected. For compound 
(3) four inter-multiplet excitations have been observed. 
All the magnetic excitations are marked in Fig. [TD1 with 
the corresponding transition energies. 

To complete our investigations of the transitions within 
the S = 12 ground-state multiplet, we additionally per- 
formed high resolution measurements of molecule (3) us- 
ing IN5 with incident wavelengths of 10.5 A (FWHM = 
13 ^eV at the elastic line) (see Fig. ITT1) . These measure- 
ments allowed us to observe transitions originating from 
the top of the anisotropy barrier. 

A further confirmation of the good assignment of the 
observed excitations is provided by the study of their Q- 
dependence. As revealed by Fig. [TJ1 the intra-multiplet 
transition (AS=0) shows a distinctive Q-dependence, 
with a pronounced intensity at low Q, that dies out 
quite rapidly following the Mn 3+ form factor. In con- 
trast, inter-multiplet excitations present flatter behavior, 
with considerably less intensity at low Q. The assign- 
ment of the observed excitations to intra-multiplet or 
inter-multiplet transitions has been confirmed by com- 
parison with FDMR measurements performed on both 
compounds (see Fig. [7] and Fig. |8]). The position of 
the intra-multiplet INS transitions are consistent with 
the FDMR measurements performed on the same sample 
(see Table HJ . Due to the different selection rules of INS 
and FDMR, we can conclude that all the peaks observed 
at T = 2 K above 1.2 meV energy transfer correspond to 
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FIG. 10: (Color online) INS spectra collected on IN5 with 
incident wavelength of 3.4 A at T = 2 K (blue circles) and 16 
K (red squares), (a) for sample (2) and (b) for sample (3) . 
The observed transitions are labeled with the corresponding 
transition energies in meV. 




-1 -0.8 -0.6 -0.4 -0.2 
Energy transfer (meV) 

FIG. 11: (Color online) High resolution INS spectra of 
molecule (3) collected on IN5 with incident wavelength 10.5 
A at 24 K. The energy gain spectra is displayed. Continuous 
lines are the calculations using the spin Hamiltonian of Eq. 
©■ 



inter-multiplet transitions, since they are absent in the 
FDMR spectra. 

The straightforward assignment of the base temperature 
observed excitations allows us to draw some consider- 
ations on the experimentally deduced energy level di- 
agram. For both compounds, a rough estimate of the 
splitting of the spin ground multiplet gives I-DIS" 2 ~ 6.5 
meV. This value is comparable to the energy interval ex- 
plored by the high energy transfer INS data (Fig. ITU1) . 
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FIG. 12: (Color online) (a) Energy-wavevector colormap of 
sample (3) collected on IN5 with incident wavelength of 5.0 
A. (b) Q-dependence of two transitions from the ground state. 
The black squares correspond to the |S = 12, M s = ±12 >— > 
\S = 12, M a = ±11 > intra-multiplet transition and the green 
circles display the \S = 12, M s = ±12 -t \S = 11, M s = 
±11 >. 



where most of the inter-multiplet S = 12 — > S = 11 
excitations have been observed. This experimental ob- 
servation leads to the conclusion that also in (2) and 
(3) several excited states lie within the anisotropy split 
ground state, with the consequent breakdown of the GSH 
approximation. Due to the inadequacy of the GSH for 

(2) and (3), the microscopic spin Hamiltonian (Eq. [5]) 
was used to model the data and extract the exchange 
constants and anisotropies . The minimal set of free pa- 
rameters is given by three different exchange constants 

J\V = Jl, J\1 = J23 = Jl3 — J\'T — JtZ 1 — Jl'3' = J2, 

and J13' = Ji'3 = J3 (Fig. \5§ and two sets of crystal- 
field (CF) parameters d\ — dy , c\ — ci<, and d-i = dy , 
c 2 = c 2 '. Indeed, the ligand cages of sites 1 and 3 are 
rather similar and we assumed the corresponding CF pa- 
rameters to be equal. Since experimental information is 
insufficient to fix independently the two small c param- 
eters,we have chosen to constrain the ratio c\jci to the 
ratio di/c?2- 

The isotropic exchange and crystal field parameters de- 
duced by the simultaneous best fit of the experimental 
data are reported in Table |TTJ Figure [13] shows the cal- 
culated energy level diagram using the best fit procedure 
for Eq. H (left) and the GSH model (right) for (2) and 

(3) . 
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TABLE II: Isotropic exchange and CF parameters for Eq. [2] 
(in meV) deduced by fitting INS and FDMR data for the two 
Mn6 S=12 compounds. 
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FIG. 13: (Color online) Calculated energy level diagram for 
molecule (2) (top) and molecule (3) (bottom). The level 
scheme on the left side is calculated using the microscopic 
spin Hamiltonian in Eq. [2] while the level diagram on the 
right has been calculated in the GSH approximation (Eq. [T|). 
The dashed lines correspond to the observed value of U e g. 



IV. DISCUSSION 

The experimental data collected on the three variants 
of Mn.6 clusters provide direct evidence that a general fea- 
ture for this class of compounds is the nesting of excited 
multiplets within the ground state multiplet. This is an 
unavoidable effect when the isotropic exchange parame- 
ters have the same order of magnitude as the single ion 
anisotropy parameters, as it happens to be for Mng. The 
nesting of spin states can be clarified by observing the 
energy level diagrams for the three molecules presented 
in Fig. [5] and Fig. [131 The diagram on the left shows 
the energy levels calculated by a diagonalization of the 
full spin Hamiltonian, while the energy level scheme on 
the right hand side has been calculated considering the 
GSH approximation. It is clear that the GSH does not ac- 
count for any of the spin states with S different from Sqs 
that lie within the split GS energy level diagram. The 
above states represent a shortcut for the relaxation of the 
magnetization and can promote resonant inter-multiplet 
tunneling processes that manifest as additional steps in 



the magnetization curve absent in the GS model£ ' 30 i 33 ~ — . 
The overall result is a lowering of the effective anisotropy 
barrier with respect to an ideal molecule where the spin 
ground state is well separated from the excited ones, as 
was firstly demonstrated in Ref. 0. 

We have calculated the relaxation dynamics of molecule 
(1) following the same procedure adopted in Ref. HI for 
molecules (2) and (3). We applied a master equations 
formalism in which the magnetoelastic (ME) coupling is 
modeled as in Ref. HH, with the quadrupole moments as- 
sociated to each triangular unit isotropically coupled to 
Debye acoustic phonons. 

The transition rates W s t between pairs of eigenlevels 
of the dimer spin Hamiltonian Eq. [3] is given by: 



W st = 7 2 A* t n(A st ) £ (s\O quq2 (S A )\t)( s \O quq2 (S B )\t) 

A,B,q u q 2 

(4) 

where O qi „ 2 (Sa,b) are the components of the Cartesian 
quadrupole tensor operator, n(A st ) = ( e hA st/k B T _ -^-l 
and A s t — (E s — E t )/h. We found out that the resulting 
relaxation spectrum at low T is characterized by a sin- 
gle dominating relaxation time whose T-dependence dis- 
plays a nearly Arrhenius behavior r = tq exp(J7/fcsT), 
as previously observed for molecules (2) and (3) 5 . The 
relaxation dynamics of M is indeed characterized by two 
separated time scales: fast processes that determine the 
equilibrium within each well of the double- well potential 
and a slow inter-well process that at low temperature de- 
termines the unbalancing of the populations of the two 
wells, and thus sets the time scale for the reversal of the 
magnetization. As can be observed from the energy level 
diagram of Fig. [5] there are several levels that can be in- 
volved in the inter-well relaxation process, giving rise to 
an overall effective barrier U c s different from the simple 
energy difference between the M = and M = ±4 states. 
The corresponding calculated energy barrier U ca \ c — 32 K 
reproduces quite well the experimental value, J7 e ff= 28 K. 
The lowering of the barrier is therefore attributed to the 
presence of these extra paths. Indeed, the calculations 
for artificially isolated S = 4 yield U = 47 K. 
It is worth commenting also on the D value for the ground 
state of each molecule. Whilst no large difference be- 
tween the local d of the low (1) and high ((2) and (3)) 
spin molecules is expected, the overall D value, as de- 
termined using the GSH approximation, is much higher 
for the S = 4 molecule (D « —0.263 meV) than for the 
high spin molecules (D « —0.045 meV ). However, this 
observation should not be misinterpreted. The difference 
arises from the fact that D depends on the projection 
of the individual single-ion anisotropics of each magnetic 
ion onto the total spin quantum number S. In the case 
where the S'-mixing is negligible and the spin ground 
state is a good quantum number, the D parameter for 
a specific state S can be written as linear combination of 
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the single-ions anisotropy tensors (Ref. [25h : 

TV 

D = ^ ai d, (5) 

»=i 

The projection coefficients of the single ion anisotropy 
to spin states of different 5 values can differ significantly, 
giving rise to considerably different D values. The lig- 
and field study of various members of the Mng family 
(Ref. l37l ) provides experimental evidence of this. Re- 
cent theoretical studies proposed that the intrinsic re- 
lationship between 5 and D causes a scaling of U that 
goes approximately with 5° (see Ref. and l26l). rais- 
ing the question whether it is worth trying to increase 
the value of spin ground state to obtain a larger energy 
barrier. Indeed, higher spin ground states would cor- 
respond to lower D parameters, neutralizing the overall 
effect on the height of the anisotropy barrier. In recently 
performed electron paramagnetic resonance studies the 
authors proposed that the barrier goes roughly with 5 1 
instead^. In the specific case of Mn6, because of the 
very large 5-mixing, the projection onto a well defined 
spin state is no more justified and it is not possible to 
associate the barrier U to a defined 5 value. However, 
if we consider the effective anisotropy barrier for artifi- 
cially isolated 5 = 4 and 5 = 12 states (i.e. U = 47 K 
for (1) and U = 105 K for (2)), we can confirm that the 
barrier does not go quadratically with 5, as one could 
naively deduce from the equation U = |-D|5 2 . Indeed, 
Us=i 2 /U s=4 = 2.2 <C 12 2 /4 2 =9. This confirms what 
has been pointed out in Ref. 0, i.e. even though the 
highest anisotropy barrier is obtained with the molecule 



with the highest spin ground state, the increase of the 
total spin is not as efficient as one would expect and al- 
ternative routes, like increasing the single ion anisotropy, 
should be considered. 



V. CONCLUSION 

We have performed INS and FDMR measurements 
on three variants of Mng molecular nanomagnets, which 
have the same magnetic core and differ by slight changes 
in the organic ligands. INS measurements have unam- 
biguously evidenced the presence of low lying excited 
states in all the three molecules. The combination of 
the two techniques enabled us to determine the spin 
Hamiltonian parameters used for the analysis of the mag- 
netic properties. The nesting of excited states within the 
ground state multiplct strongly influences the relaxation 
behavior and plays a crucial role in lowering the effective 
energy barrier. The calculations of the relaxation dynam- 
ics give results that are consistent with the experimental 
values and show that the highest barrier is obtained for 
ideal molecules with an isolated ground state. This ob- 
servation might be valid for a wider class of SMMs. 
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